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I. INTRODUCTION 



Spintronics-based quantum computation systems are expected to be one of the important 
successors of the microelectronic-based conventional computation systems in the future. 
The essential processes of realizing spin computation, spin injection into devices, and spin- 
related transport in semiconductors have attracted enthusiastic attention in the past few 
years. To avoid the practical difficulty of integrating devices, some electrical methods, 
instead of conventional ones applying an external magnetic field or employing circularly 
polarized light, are required to induce spin-polarized carriers in semiconductor devices. One 
simple idea is to use the ferromagnet-semiconductor interface to produce spin-polarized 
electrons, but this method must face the mismatch of physical parameters between these 
two quite different materials |T|J5|. The employment of diluted magnetic semiconductors 
(DMS), which can match well with other extensively used semiconductors like AlGaAs, have 
provoked a lot of interest in DMS Recently, based on the Rashba spin-orbit interaction, 

an intrinsic effect in inversely asymmetric or asymmetrically confined nanostructures of 
non- magnetic semiconductors (NMS), several designs have been proposed to spin-polarize 
electronic currents in nanostructures This progress in spintronics offers the possibility 
of doing the spin injection in conventional materials and is bringing more and more focus 
on how to control and utilize the Rashba effect in these well-known materials and well- 
controlled structures. The possibility of establishing a spin transistor, based on the Rashba 
interaction, has also been considered, but further investigation is required to obtain devices 
of good behavior. 

Spin degeneracy of carriers in semiconductors is a result of inversion symmetry, in space 
and time, of the considered system. By introducing a spatial inverse asymmetry, one can re- 
alize spin splitting for carriers of finite momentum, without applying any external magnetic 
field. This so-called Rashba spin-orbit interaction P,|T0f has been confirmed experimentally 
in different semiconductor structures [pTH I^|. In semiconductor heterostructures, this spa- 
tial inverse asymmetry can be easily obtained by either built-in and external electric fields or 
by the position-dependent band edges. It is found that in many cases, especially in narrow 
gap semiconductor structures, the corresponding spin-orbit interaction is a linear function 
of the electronic momentum k expressed as the Rashba term !? ■ (k x E) in the electron 
Hamiltonian, where !? is the Pauli spin matrix and E the local electric field. Thus, a local 
electric field works on the electronic spin like a local magnetic field perpendicular to the 
directions of the electric field and of the electron momentum. The averaged Rashba param- 
eter is proportional to the average electric field weighted by the electron probability and 
can be well controlled by a top (back) gate over (below) the device. Recently, Nitta et al. 
I7)| , [TE|] studied the dependence of the spin-orbit interaction on the surface electric field in 



an inverted InGaAs/InAlAs heterostructure and Grundler (T^j showed that the penetration 
of the electron wave function into the barrier layer can greatly enhance the spin-orbit inter- 
action in InAs quantum wells because in this case the electrons suffer a stronger effective 
electric field due to the band-edge difference. The Rashba effect began to be considered 
as one of the powerful tools in making spintronic devices after the pioneer proposal of the 
spin-polarized field effect transistor by Datta and Das [18|]. Mireles and Kirczenow [fl9 



studied in detail the ballistic spin-polarized transport and the Rashba spin precession in a 
semiconductor waveguide using a tight-binding model but used somewhat large values for 
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the strength of the Rashba parameter. Based on the fact that the spin orientation of elec- 
trons in nanostructures depends on the direction of their momentum, Kiselev and Kim 
proposed a T-shaped spin filter while Governale et al. pi announced recently the possibility 
of making a more effective spin filter with tunnel-coupled electron waveguides. Spin injec- 
tion into non-magnetic semiconductors with the help of magnetic metals [0,^O| or diluted 
magnetic semiconductors |J has also been intensively studied. However, it is not yet known 
how to effectively control the spin-polarized flux in those ballistic transport devices and 
waveguides, though ballistic electronic transport, disregarding spin polarization, has been 
studied in detail in the past several years pi| - |25|] . In a previous paper we showed briefly 
how spin-polarized transport can be produced and controlled in stubbed waveguides when 
only spin-up or spin-down electrons are injected and only one mode is propagating in the 
waveguide. 

In the present paper we study in detail the possibility of a spin transistor which can 
control the flux strength (transmission rate) and spin orientation using periodically stubbed 
semiconductor waveguides in which the Rashba effect is present. In addition, we consider the 
case when two modes are allowed to propagate in the waveguide as well as that of " defects" 
in a periodic array of stubbed waveguides. We further consider the case of injected electrons 
polarized partially up and partially down and new stub shapes as well as the influence of 
finite temperatures on the conductance. The results obtained are mentioned in the abstract 
and are detailed as follows. In Sec. II we present the formalism and in Sec. Ill the numerical 
results for one (III. A) and two propagating modes (III.B) with spin-up injection. Results 
for finite superlattices with defects are considered in Sec. III.C and the injection of electrons 
with their spins polarized in an arbitrary direction in Sec. III.D. Conclusions follow in Sec. 
IV. 



II. FORMALISM 



For a typical two-dimensional (2D) electronic system in the x — y plane in narrow 
gap semiconductor nanostructures such as InGaAs/InAlAs quantum wells, the one-electron 
Hamiltonian including the lowest order of spin-orbit interaction can be expressed as 



H 



2D 



+ a ( " x ~) 
2m* h 
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+ ia(a y — 
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(1) 



where V 2 = d 2 /dx 2 + d 2 /dy 2 , and V ± = d/dx ± id/dy. The parameter a measures the 
strength of the spin-orbit coupling and is proportional to the interface electric field; a = 
(a x , <J y ,<J z ) denotes the spin Pauli matrices, and p is the momentum operator. In the presence 
of a we assume that the new wave function has the form 



${k x , k y ) = e ikxX+ik y y C a \a) = e lkx 



x+ik y y 



c- 



(2) 



with | a) = ( J ) (spin up) or ( \ j (spin down). The solutions of the equation H 2D ^(k XJ k y 
E^/(k x , k y ) are readily obtained as 
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the corresponding eigenvalues are 

E {k X i ky 



1 

V2 



1 



2m* 



k 2 ± ak, 



(3) 



where k = ^Jk^ + k%. The electrons are now spin polarized and oriented perpendicular to 
the electronic momentum in the 2D plane. 

If the electron gas is confined along the x direction by a potential V(x), such as the 
one in the stubbed waveguide shown in Fig. [TJ(a) , we have a quasi-one-dimensional (Q1D) 
electronic system. The Hamiltonian becomes 



H 1 V72 

-aV 



+ V(x) 



-&^ 2 + V{x) 



(4) 



Denoting by (f) n (x) the solutions of the equation — (h 2 /2m*)d 2 /dx 2 + V(x) (j) n (x) 
E n (p n (x), we can express the Q1D eigenstates in the form 



C: 



(5) 



in each of the regions I, II, or III we have <f> n {x) = sm(mr(x + w/2)/w), where w is the width 
of the region along x. Then the equation ifQ 1D \|/ = E^ takes the form 



E 



[E n + h 2 k 2 y /2m* - E\<\> n Ct + ak x (j> n C- + a^C' 



^ \ akJ n C+ - a<jJ n C+ + [E n + h 2 k 2 /2m* - £]0 n C~ 
where <p' n = d<fi/dx. 



0. 



(6) 
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FIG. 1. (a) Schematics of a stub tuner with two units. The waveguide segment between the 
two stubs (shaded regions of width b) has length I and widths a and c as indicated. The asymmetry 
parameter d is the distance between the lines, along the y axis, that bisect the main waveguide and 
the stubs as indicated, (b) Dispersion relation for a waveguide based on Eq. (^). The horizontal 
dashed line indicates the four possible values of the wavector for the same energy E. 



Mutiplying both sides with m (x) and integrating over x leads to (/ dx<f) m (x)4> n (x) = 5, 
( [E m + h 2 k 2 /2m* -E}C+ + ak x C~ +E n f dxa^' n Cn " 



ak x C+ -En I dxa<p m <P' n C+ + [E m + h 2 k 2 /2m* - E]C~ 



0. 



(7) 
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According to degenerate perturbation theory, if the inequality 

< 1 

\. Then Eq. ([7p becomes 

= (8) 



holds, we can neglect the subband mixing term / dx<ft r , 



E m + (h 2 /2m*)k 2 - E 



ah, 



E m + (h z /2m*)k 2 - E 




and its eigenvalues are 



E ± (k y ) =E m + (h 2 /2m*)k 2 ± ak y 



(9) 



The eigenvectors corresponding to E + ,E~ satisfy = ±C~. Accordingly, the spin eigen- 
functions are taken as |±) = ±i )• The dispersion relation E ± (k y ) versus k y resulting 
from Eq. @ is shown in Fig. 0(b). Note that the minimum of subbands are shifted from 
k y = by the value A/2 = m*a/h 2 . In each subband, electrons of the same energy E have 
four different momentum values k y , i.e., the positive or negative ky and k~ values pertain- 
ing to the branches E + (k y ) and E~{k y ). In fact, electrons belonging to different branches 
have opposite spin orientation in the x direction. The difference in wave vectors ky and k~ , 
resulting from E + = E~ = E, reads 



k y k y 



2m* a /h 2 = A. 



(10) 
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Now let us consider the transmission process when an electron of energy E is incident 
from the left to a stubbed waveguide as illustrated in Fig. |l](a). The electron wavefunction is 
decomposed into the plus ( } ) and minus ( _\ ) branches and the procedure outlined above 
applies to each of the three regions labelled I, II, and III in Fig. [l|(a). In each region we 
have 4> n (x) = sin(n7r(x + w/2)/w), where w is the width of the region along x. Including 
spin and referring to Fig. |l|(b) we can write the eigenfunction of energy E in region I as 



E J + <W 



-1 



and in region III as 



x sin(c m (x + c/2)) (11) 



x sin(a n (x + a/2)); (12) 

here c m = rrm/c, (3 m = (2m*E — c^) 1//2 , a n = nn/a, and a n = (2m*E — a 2 /) 1 / 2 where 
m, n — 1, 2, 3, ■ • • denote the order of the transverse modes. The symbols a" m , b\ m , a% n , b% n 
(a = ±) represent the coefficients of different electronic modes existing in the device. Similar 



to the procedure of matching electronic wavefunction when spin is disregarded [22], in region 



II we use two auxiliary sets of solutions to the wave equation, one of which matches the wire 
on the left and the other on the right, with each vanishing elsewhere on the boundary: 
0s — J2ki.Xk + Xk)- The appropriate boundary conditions are 

x *(y = 0,x) = 0; 

0, x < -a/2; 

\!:(!l = b..c)= { sm[a k (x + a/2))("f + y + ?| + J 2 > V -a/2 < x < a/2; (13) 



a 2k + a 2k b 2k b 2k 



x > a/2; 



0, x < -c/2; 



XUV = 0,x) = sin[c fe (x + c/2)] \ k ' l k A* u- , -c/2 < x < c/2; (14) 

V a ik ' a ik °ik °ik ) 

0, x > c/2; 



and Xk (v = b, x) = 0. Xk can be expanded as 

x R k = E ( J ) + v^^y ' 1 



l 



u fl e -i(TH+A)y })_ v R e -i7nV ( 1 



"n ^ \ 1 / \ 1 

x sin(/i n (x + /i/2 - d)). (15) 
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Multiplying Eq. fll3|) by sin(a m (x + h/2 — d)), integrating over x (—a/2 < x < a/2) and 
using Eq. (|l5l), we obtain 



(16) 



where 



/•a/2 

-^fcL — / ^ sin[afc(x + a/2)] sin(a m (x + h/2 — d), if = h[e 

J -a/2 



and 7 m = (2m* E — a^) 1 / 2 . Similarly, x£ can be expanded as 



1 



^ e -iTn(2/-6) 



where 



and 
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M m = ^(a^ + bf k )lj/ m / K m , t> m = — 2(a lk + b lk )I km J 'K+ 



km 



c/2 
-c/2 



dx sin[cfc(x + c/2)] sin[a m (x + a/2 — d)] 



(18) 
(19) 

(20) 



Requiring the continuity of the derivative of the wavefunction at the interfaces, <pi and 
2 must satisfy the following equation: 
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(21) 



The elements of sub-matrices A, B, C, D, E, F, G, H, a, (3, a', /3' have values 
A lk = E4(2 7m + A)l£ n I* a /{cK+), 

m 

B lk = E4[7 m e"^ 6 + ( 7m + A)e^ +A)b }lLlL/(cK-), 

m 

C lk = E -4(2 7m + A)ltlL/(cK-), 

m 

E "4[7 m e^ b + ( 7m + ^)e~ i ^ +A) VLlL/(cK+), 

m 

E4[7 m e^ b + ( 7m + A) e -^ +A ) 6 ]/«/ fe B m /(air+), 

m 

E4(2 7m + A)/^4 L m /(aK-), 
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-i-/mb 



E-4(2 7m + A)/^4 L m /(air m ), 



(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 
(29) 



where = /3 k 8i k , a ik = a k 5 ik , j3' lk = {j3 k + A)5 ik , and a' lk = (a k + A)5 ik . 

Together with the matrix P corresponding to a waveguide segment of length I rather 
than to a stub, 
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(30) 



we can connect the incident waves (to the left of region I) with the outgoing ones (to the 
right of region III) via a transfer matrix M 



( at \ 

bfn 
V Kn ) 
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btut 



a. 



out 

V Kut J 



(31) 



Here M is a 4 x 4 matrix with sub-matrices M mn as its elements. If we assume that a 
spin-coherent electron beam of single energy is injected into the device from the left and 
it is detected at its right end by an analyzer, there is no backward propagating wave at 
the output and b^ ut = b~ ut = 0. The coefficients of spin-up and spin-down electrons are 
af n + a^ n and af n — a~ n for the input wave and 2+^ + a~ ut and a^ ut — a~ ut for the output wave, 
respectively, if we neglect the injection mismatch of electronic momentum between different 
electron branches. Then the spin-dependent transmission rate can be calculated similar to 
the procedure used when the spins are disregarded. 

In the design of a spin transistor first proposed by Datta and Das, we can introduce 
the above stubbed waveguide together with the Rashba effect into the device to control and 
flip the spin current. To see this effect clearly, we connect a spin polarizer (analyzer) to 
the left (right) of the structure and inject spin-up polarized electrons into it and detect the 
polarization of the outgoing electrons. The electron beam is equally decomposed into the plus 
branch { ) and the minus branch ^ \ ^ if the momentum difference between two branches 
can be neglected when considering the wavefunction match of input and output. Since the 
electrons belong to different branches will get different phases when propagating along the 
device, the output spin orientation of the electrons, which is a result of the composition 
of electronic wavefunction of the two branches when leaving from the waveguide, can be 
different from the initial one. In our case, the spin-up transmission rate T + and spin-down 
transmission rate T~ read 
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(32) 



where a out - m n denotes the mth mode iuput contribution to the nth mode output a outt n and 
the sums over m,n is over all possible propagating modes (number N m ) in the waveg- 
uide. For spin-up electrons injected into a simple waveguide of length /, we can as- 



sume a; 



l out,mn 



in,m 



1/2 and the output coefficients can be calculated using Eq. (^ 



and a 



out,mn 



e ^ am+A ^ l 8 mn . Since a n = (3 n for a simple waveguide, we 



S 



have = N m 1 ± e lAl . The probability of detecting a spin-down ( ° ) electron will be 

T~ = N m sm 2 (S9/2) and that for a spin-up ( J ) electron T+ = N m cos 2 (56»/2), where 

50 = A/ is the phase difference between the two spin propagation modes |+) and |— ) in 
the structure after traveling a distance /. In stubbed waveguides, it will be shown that the 
percentage of transmitted spin-up and spin-down electrons follows that same rule as in a 
stubless waveguide when the Rashba effect is weak. The advantage of stubbed waveguides is 
that using side gates to control the length h of the stubs and the distance d of their centers 
from that of the waveguide, we can control the transmission rate. At the same time we 
can adjust the back gate bias to change the Rashba parameter and then the outgoing spin 
orientation. 

Once is known, the conductance G = G(E, 0) at zero temperature is given by 
the Landauer-Biittiker formula, G ± = (2e 2 /h)T ± . For finite temperatures the conductance 
G(E, T) is given by 



G(E, T) = jT G(e, T = 0) (-f) de, (33) 



where f(e — E) is the Fermi function. 

In real devices the shape of the stubs can be different than a rectangular one and one 
question we need to answer is the influence of the shape of the stubs on the transmission 
output. It is known from electronic stub tuners that stubs of different shape, e.g., Lorentzian 
or triangular stubs, do not change the qualitative behavior of the transmission. To quantita- 
tively study the transmission rate of electronic current through a waveguide with arbitrarily 
shaped stub, we can break the stubs into a series of rectangular segments with the same 
width hi and different heights. Each segment is described by a transfer matrix Mj. The 
complete shape is well described by the product M T = flj Mi if the segment width is much 
smaller than the electronic wavelength. The entire procedure described above for a single 
stub is then repeated each segment as many times as required by the particular shape. 



III. RESULTS AND DISCUSSION 

In the numerical calculations we consider a stubbed Jno.53Gao.47As/JnA/As waveguide. 
The effective electron mass is m* = 0.042mo and a typical average spin-orbit constant 
a = 1.6 x 10 _11 eVm is assumed throughout the paper unless otherwise specified. To verify 
the validity of the perturbation theory we evaluated the bound states of one unit, of stub 
length h = 2000A and width b = 150A connected to a waveguide of width a = 250A and 
segment length I = 100A to the left and right of the stub. We find that the ratio of the 
intersubband mixing energy over the difference between the lowest two bound states is less 
than 10%. We also calculated the energy bands of a superlattice of such stub units with a 
waveguide segment of length I = 200A between two consecutive stubs. The separation of 
the two spin bands is less than 15% of the band energy. 



A. One mode allowed in the waveguide 
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1. Rectangular stubs 



To view the parameter dependence of the spin transmission through a stubbed waveguide, 
we inject spin-up polarized electrons into a waveguide of width a = 250A with a symmetric 
double stub of width b = 150A and measure the output flux rate of spin- up and spin-down 
electrons. Considering only the first mode and for a waveguide segment of length I = 450A, 
we plot in Fig. || the spin-up electron transmission as a function of the stub length, reduced 
by the width of the waveguide (h — a), and of the electron energy E. Only electrons 
with energy higher that the first subband in the waveguide Ei = it 2 / (2m*a 2 ) = 14.3meV 
can pass through the device and the output percentage of spin-up and spin-down electron 
depends only on the total length of the device. Here A = 2m*a/h 2 = 0.1764 x 10 8 m _1 
so the maximum of the spin-up electron transmission rate is cos 2 A(b + I) = 76.4% and 
the maximum of the spin-down electron transmission rate 23.6%. The stub begins to play 
a role in adjusting the transmission rate when electrons have an energy higher than the 
first subband in the stub E{ = n 2 /(2m*b 2 ) = 57.2meV. If the electronic energy has a 
value between the first and second transverse subbands, i.e., if the relation E[ < E < 
E% = 2-7T 2 / (m*b 2 ) holds, only one transverse mode enters the stub and a simple transmission 
pattern appears with transmission gaps along the curves h = n\ + h , where ho is the 
position of the first gap and n = 0, 1, 2, • • •. Notice that the period A is different from the re- 
direction wavelength of the first subband wavefunction in the stub Ai = 2rcb / \f2m* Eb 2 — tt 2 
and depends on the parameters a and I. When the width b of the stub becomes wider, more 
modes can exist in the stub and the transmission pattern becomes complicated, see Sec. 
IIIc. 




FIG. 2. Three-dimensional illustration of transmission T versus electron energy and stub arm 
length h — a of a symmetric stub tuner. 
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FIG. 3. Spin-down transmission of an array of five stub tuners as a function of the stub length 
h and stub shift d when only spin-up electrons are injected. Cross sections along the dashed lines 
are shown in Fig. (|j). 

We can also change the asymmetry parameter d, the distance along x between the mid- 
point of the waveguide and that of the height h shown in Fig. 1(a), to shift the stub along 
the x direction so the structure becomes asymmetric. In Fig. |3| we show by a wired surface 
the transmission rate for electrons of fixed energy E = 48meV as function of d and h for a 
five-stub device with parameters a = c = 250A, b = 150A and I = 207. 5A. The parameter 
I is chosen to satisfy cos 2 [5A(6 + /)] =0 so that spin-up electrons are totally blocked, ( 
T + = 0), and only spin-down electrons come out. The transmission is a periodic function of 
d and h and the gaps appear in the triangle-shaped regions. The centers of these triangles 
are located at the points in the h — d plane that satisfy h = nX/2 + ho and 2d + h = nX + h 
for integer n — 0, 1, 2, • • •. In this figure we find A = 558 A, ho = 694A and Ai = 660. 5A. 
One of the interesting facts found in Fig. |3] is that a square type transmission curve with a 
wide gap can be obtained by changing d for a fixed value of h. 
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FIG. 4. Two-dimensional cross section of Fig. (]3|) in the transmission T vs stub length h for 
fixed d ( d = in a) and d = A/8 in b) ) and vs stub shift d for fixed h ( h=1289.5A in c) and 
h=1568.5A in d) ). 

To clearly show when and how these square-shaped transmission curves appear, some cross 
sections of Fig. [| are shown in Fig. |4j. In Fig. |4](a) the spin-down transmission rate is shown 
as a function of h for the symmetric case (d = 0). Gaps appear with period A and the first 
gap begins at h — ho — 694A and has a width 5h = 4fjA. If we cut the surface along h at 
d = 69.75A = A/8, we obtain Fig. |(b). Here we see a periodic structure similar to that in 
(a) but with half the period A/2 and half the gap width 5h = 20A. If we fix the value of 
the stub length h in one of the gaps in Fig. ^(a) and cut the surface of Fig. |3] along d, we 
obtain the curve shown in (c), where h = 1289. 5A ~ h + A. The transmission has a period 
A/2 along d with gaps as wide as 8d = 91 A. The same periodic transmission versus d curves 
is shown (d) when we make the same cut as above but fix h = 1568. 5A ~ h + 3A/2 so the 
gaps appear at d = A/4 rather than at d = 0. We see that in (c) and (d) we can completely 
flip or block the input electronic spin flux and the transmission has an almost square-wave 
dependence on the adjustable parameter d with the gap/band ratio as big as 0.5. 



2. Stubs of general shape 

To study the influence of the stub shape on the transmission output we proceed as 
outlined at the end of Sec. II. As an example we consider devices with double stubs of 
Gaussian shape. As a function of x the y coordinate of the boundary of each stub is taken 
as y — (h — a)/2 ± d](e~ 8x ' b — e~ 2 )/(l — e~ 2 ). In Fig. |^(a) we show the transmission of 
a structure of two symmetric double stubs as a function of the total length h of the stubs. 
Only spin-up electrons are injected. The thin and thick solid curves show, respectively, the 
spin-up and spin-down transmission for Gaussian-shaped stubs and parameters a = 250A, 
b = 375A, / = 275A; the thin and thick dotted curves are for rectangular with parameters 
a = 250A, b = 290A, L = 360A. We see that the transmission through the waveguide with 
Guassian-shaped stubs is similar to that with narrower rectangular stubs for the same total 
length (b + L). Such a similarity was noticed earlier between rectangular and triangular 
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stubs . If we fix the total length of the stubs of both waveguides h = 825A and shift the 
stubs, we obtain Fig. b)) for the transmission as a function of the shift d. 




h(100A) d(100A) 



FIG. 5. Transmission as a function of the stub length h, for symmetric double stubs (a), and of 
the stub shift d for asymmetric double stubs (b). The thin and thick solid curves show, respectively, 
the spin-up and spin-down transmission for Gaussian shaped stubs. The thin and thick dotted 
curves show the corresponding results for rectangular stubs with parameters a = 250A, b = 290A, 
L = 360A. 



B. Two modes allowed in the waveguide 

When the electron energy is high enough, two or more transverse modes can propagate in 
the waveguide. In Fig. |6|(a), we plot the transmission as a function of the electronic energy 
E through a device with one waveguide segment of width a = 250A and length / = 250A and 
one stub of width b = 125A and length h = 1250A. Here the total device length is adjusted 
so that both spin-up (solid line) and spin-down (dotted line) electrons can be observed at 
the right end. When the electron energy is lower than E% = 57.2meV, only one propagating 
mode exists and the transmission is almost flat; the total transmission (T + + T~) approaches 
unity because in this energy region the stub is too narrow to give significant contribution. 
With the increase of the electron energy, two modes can exist in the structure. In some 
cases, two modes compete and a deep minimum exists with the transmission dropping from 
2 to a value much lower than 1 as shown in Fig. 0(a). If a wider stub is used here, more 
modes can exist in the stub and the curve has many more oscillations. 
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FIG. 6. (a) Spin-up T + (solid line) and spin-down T~ (dotted line) transmission through one 
stub tuner as a function of the electron energy E. (b) T + and T~ for four identical stub tuners as 
a function of the stub length h. 

In Fig. ||(b) we inject electrons of energy lOOmeV into a device of four equal units, each 
of which has one asymmetric stub with b = 125A and d = 30A connected to one waveguide 
segment of width a = 250A and length I = 300. The transmission rate is shown as a 
function of the length h of the stubs. In this figure the spin-orbit parameter is adjusted to 
a = 1.8 x 10 _11 eVm so that almost all electrons are spin-flipped after passing the device. As 
in Fig. 4 we observe a series of gaps, narrower than the corresponding ones of Fig. 4, but 
now with two modes present and electrons of either spin completely blocked. Notice that 
the transmission rate is still a periodic function of the stub length h despite the presence of 
two modes. 
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C. Introducing defects 




FIG. 7. Spin-up transmission T + versus electron energy E for a periodic array of seven units 
(solid line) and for the same structure with the central (fourth) unit replaced by a defect unit. 

Using a series of periodically arranged stub units, we can make a spintronic crystal 
(superlattice) with a transmission similar to that shown in Fig. |3|, where five units are used. 
In this kind of crystals, we can expect to obtain devices of better behavior than that of only 
one unit. If we add some impurities into the crystal, we can change the properties of this 
crystal and observe a different transmission rate. The solid line in Fig. [7] shows the first mode 
spin-up transmission rate versus electron energy through a pure finite crystal of seven equal 
units with symmetric stubs and parameters a = 250 A, b = 250A, I = 125A, and h = 2250A. 
The transmission shows gaps near E = lOmeV, 40meV and 60meV. When we substitute 
the middle (fourth) unit by a defect, that is, a unit which has stub length hd = 4500Aand 
otherwise the same parameters, we find that the defect introduces resonant peaks in the 
gaps at E = 40 meV and E = 46 meV and E = 55 meV. The spin-down transmission rate 
has the same structure but different scale. A similar behavior was reported for photonic 
tuners in Ref. 



D. Injection of spin-up and spin-down electrons 

Up to now we considered electronic conductance or transmission through stubbed waveg- 
uides when only spin-up electrons are injected at zero temperature (G(E,T = 0)). In this 
subsection we discuss what happens if arbitrarily spin-polarized electrons are injected into a 
stubbed waveguide at finite temperature. In the a z spin representation the spin-orientation 
is denoted by the coefficients C + and C~ . For instance, ( ±i ) describes a spin oriented 

along the x-direction whereas f ±» ) describes a spin oriented along the y-direction. In Fig. 
^|(a) and (b) we show the spin-up (dotted line) and spin-down (solid line) conductance for a 
structure composed of eight stub units with a = 150A, b = 80A, I = 69A and Fermi energy 
Ep = 133meV. The spin orientation of the incident electrons is as follows: we let C + = 0.7 
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and C = —0.3; since the state must be normalized (|C + | 2 + \C | 2 = 1), this corresponds 
to 84.5% spin-up and 15.5% spin-down incident electrons. 




h (100A) d (100A) 



FIG. 8. Spin-up G + (lower lines) and spin-down conductance G~ (upper lines) as functions of 
the stub length h (a) and of the stub shift d (b). The various curves and temperatures are specified 
in the text. 



As a function of the stub length h the conductances G + and G~ , at T =0 and 8K, are 
shown in Fig. H(a); the solid and dotted curves are for G + and the thin and thick solid lines 
for G~~ . As a function of the shift d G + and G~ are shown in Fig. §(b) at T =0, 4.2, and 
8K. The solid, dashed, and dotted curves are for G + and the thin, medium, and thick solid 
lines for G~ . Here we have 15.5% of the output electrons having spin up and 84.5% of them 
spin down. The conductance G~ shows a behavior similar to that of the transmission in 
Fig. |](a) and (d). As can be seen, finite temperatures smoothen the curves obtained at zero 
temperature similar to the case of spinless electrons [p2| |. As the ratio Ep/ksT increases the 
curves are smoothened or rounded off more strongly. 



IV. CONCLUDING REMARKS 



Employing the transfer-matrix method we combined the spin precession in a waveguide, 
due to the spin-orbit coupling, with the basic physics of a stub tuner, and applied it to 
the ballistic spin transmission through periodically stubbed waveguides. We found that the 
spin polarization can be well controlled by adjusting the total length of the device for a 
wide range of the electronic energy. In particular, we showed that, given spin-polarized 
electrons injected into a stub structure, we can select the spin of the outgoing electrons 
to be the same as or opposite to that of the injected electrons. We demonstrated this for 
stubs of rectangular or Gaussian shape but also for triangular stubs |Z?J. In general these 
results hold for stubs of any shape. The latter does not affect the qualitative behavior of 
the transmission but only its period when several stubs are combined. 
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More important, we saw clearly that, as a function of the stub height h and the asymme- 
try parameter d, we can have a nearly binary square-wave transmission (spin-valve effect) 
for either spin orientation, with wide gaps for stubs of different shapes and a well-controlled 
range of the stub parameters. In this respect asymmetric stubs give the best results. Their 
asymmetry, i. e., h and d, can be controlled by lateral gates in principle, such gates 
should allow for a more detailed control of the overall stub shape. A further modulation 
can be achieved if we combine several groups of stubs with different values of the spin-orbit 
coupling strength. These findings should facilitate the experimental realization of the spin 
transistor. 

A qualitative understanding of all these results is as follows. The spins precess in a 
single waveguide [Tj| due to the spin-orbit coupling. On the other hand, in a stub tuner 
waves reflected from the walls of the stub, where the wave function vanishes, may interfere 
constructively or destructively with those propagating in the main waveguide and result, 
respectively, in an increase or decrease of the transmission |[24|| . Refining this idea, it was 
shown in Refs. [[H]]- |22| that using double stubs the transmission of spinless electrons could 
be blocked completely using asymmetric stubs. Combining several stubs would result in a 
nearly square- wave transmission output, especially as a function of the asymmetry parameter 
d. The transmission shown, e.g. in Figs. is simply the result of this behavior when 
combined with the spin precession due to the spin-orbit coupling since the length of the 
device was chosen such that spin flip would occur in the stubless waveguide. 

The most clear results or simplest transmission patterns are obtained when only one 
mode is allowed to propagate in the main waveguide. If more modes can propagate in the 
main waveguide, generally the transmission pattern becomes more complex or even irregular. 
However, as we demonstrated, we can have a simple periodic transmission pattern even when 
two modes are allowed. This occurs when the stub width is short enough that only one mode 
enters in the stub region. 

Further, we showed that the above results hold when the injected spins are polarized in 
an arbitrary direction, partially 'up' and partially 'down'. We also showed that introducing 
'defects' in a finite superlattice leads to a further modulation of the spin current since 
the 'defects' produce new transmission resonances or antiresonances. One could use such 
'defects' to achieve a specific control of the transmission. 

Finally, we have seen that the effect of finite temperatures is to smoothen the zero- 
temperature conductance as in the case of spinless electrons [|^] . The degree of smoothness 
depends mainly on the ratio Ep/ksT. 
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